CN| 



O: 

< 



(N 
> 



c^ 



GLOBAL ATTRACTOR AND ASYMPTOTIC DYNAMICS IN THE 
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f^ , Abstract. We study the dynamics of the large A'^ hiriit of the Kuramoto model of 

^\] ■ coupled phase oscillators, subject to white noise. We introduce the notion of shadow 

inertial manifold and we prove their existence for this model, supporting the fact that 
O ■ the long term dynamics of this model is finite dimensional. Following this, we prove that 

.^^ ■ the global attractor of this model takes one of two forms. When coupling strength is 

\ below a critical value, the global attractor is a single equilibrium point corresponding 

C^ . to an incoherent state. Otherwise, when coupling strength is beyond this critical value, 

the global attractor is a two-dimensional disk composed of radial trajectories connecting 
a saddle-point equilibrium (the incoherent state) to an invariant closed curve of locally 
stable equilibria (partially synchronized state). Our analysis hinges, on the one hand, 
upon sharp existence and uniqueness results and their consequence for the existence of 
pH . a global attractor, and, on the other hand, on the study of the dynamics in the vicin- 

ity of the incoherent and coherent (or synchronized) equilibria. We prove in particular 
non-linear stability of each synchronized equilibrium, and normal hyperbolicity of the 
set of such equilibria. We explore mathematically and numerically several properties of 
the global attractor, in particular we discuss the limit of this attractor as noise intensity 
decreases to zero. 
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H . 1- Introduction 



The Kuramoto model is a classical model of interacting phase oscillators, coupled through 
a mean-field term, widely used to investigate synchronization. It has been applied in var- 
ious fields such as physics, chemistry and biology (see review [Ij and reference therein). 
In different contexts the same model is also referred to as the Sakaguchi model |25| . the 
Smoluchowski model [3 El [27], or the mean-field classical spin XY model l7^. One re- 
markable feature of this model is that, for large populations, it undergoes a transition 
from an incoherent state to a (partially) synchronous one as the coupling strength between 
oscillators is increased. Many numerical and theoretical studies have analyzed this tran- 
sition under various hypotheses [25]. In the present work, we are interested in the long 
term dynamics of infinitely many interacting identical noisy phase oscillators. For such an 
infinite population, the following mean field model can be derived from the microscopic 
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description (see |3| and references therein) 

' dtq{t, 9) = \dlq{t,e)- Kde [q{t,e){J * q){t,e)] t > 0, e [-vr, ^] , 

l{t,-^) = q{t,n) t>0, 

^deq{t,-Ti) = deq{t,'K) t > 0. 

/TT 
s\n{ip-e)q{t,ip)dLp, (1.2) 

i^ is a real constant representing the coupling strength (7 > such that Lq{t,9)d6 = 1. 
The variable 6 £ [— vr, vr] accounts for the phase of the oscillators, and the unknown q{t, 9) 
for the density of oscillators at phase 9 at time t. We have assumed that oscillators are 
homogeneous, i.e. they all have the same intrinsic frequency Ci; = 0, and without loss of 
generality we have assumed that the intensity of the noise perturbing each oscillator is one. 

For this system one encounters two distinct regimes depending on the ratio between 
the noise intensity and the coupling strength. When the coupling strength K is smaller 
than a critical value Kc the noise dominates, a uniform state is the only equilibrium of 
(jl.ip . and the population always tends to this incoherent state. When K > Kc instead 
the coupling dominates, a family of non-trivial coherent (or synchronized) equilibria exists, 
and the population tends to synchronize. When K > K^ for a different coupling potential 
J * q{t,9) = J^ sin'^{ip — 9)q{9)d9, existence and uniqueness of a one dimensional circle 
C = {f*{- + ip), 99 € [0,2vr[}, where/* is a non uniform probability on [0,2vr], of non-trivial 
equilibria has been established in [TS]. In a more general setting, estimates on the number 
of equilibria, and asymptotics in the large coupling limit K — )• -|-oo, have been established 
in [6]. For equation (II. ip with (II. 2p . when K > Kc, existence and uniqueness of a circle of 
non-trivial equilibria C has been established in [8j (see also |[4] and references therein). A 
spectral gap estimate and linear stability of trivial equilibria q{- + <f) have also been shown 
in [4]. 

Several major advances have recently been made in the understanding of the global dy- 
namics of Kuramoto models, in the case of deterministic oscillators and in the case of noisy 
oscillators. For an infinite number of identical oscillators subject to no noise, equivalent 
to equation (jl.ip with no diffusion, Ott and Antonsen have shown thanks to a well chosen 
ansatz that a two-dimensional invariant manifold exists, and they have provided an explicit 
expression for the orbits on this manifold |21] . These orbits are the paths followed by the 
oscillator population during synchronization, as they describe the trajectories that connect 
the incoherent state to synchronized ones. Generalizing these ideas to finite dimensional 
systems, Mirollo and Strogatz [20j have shown that regardless of the (finite) number of 
oscillators, the dynamics of the homogeneous Kuramoto model can be reduced to three di- 
mension, through an action of a three-dimensional Moebius group on the torus T^. While 
these constitute important breakthroughs in our understanding of the dynamics of coupled 
phase oscillators, they do not deal with the situation where noise perturbs the dynamics 
of the units. 

In the case of an infinite number of noisy identical oscillators (jl.ip . Vukadinovic has 
established in [27] the existence of finite dimensional invariant exponential attractors (in- 
ertial manifolds) in the invariant subspace of symmetric solutions q{t, 9) = q{t, —9) of (II. ip . 
Inertial manifolds existence theorems do not apply directly to equations of this form, but 
Vukadinovic has developed methods to show their existence for a Smoluchowski equation 
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on the circle |27| and on the sphere |28) , and for a Burgers equation [29] for example. 
Although the dynamics of (jl.ip is infinite dimensional, inertial manifolds show that it has 
some typical properties of finite dimensional systems. 

In this work we are interested in infinite populations of noisy coupled oscillators (jl.ip and 
the finite-dimensional behavior of its dynamics. This work is in line with previous findings 
and it extends them to give a complete rigorous description of the long term dynamics 
of the Kuramoto model (jl.ip . This paper is organized so as to move progressively from 
general results concerning the solutions of equation (II. ip to more refined descriptions of the 
dynamics on its attractor and ending with a number of numerical investigations opening 
the way for some conjectures. 

In section [21 we start with sharp existence and regularity results for solutions of (II. ip . 
Existence and uniqueness of solutions for (11. ip in Lp' or Sobolev H^ spaces is a classical 
result. Combining a result of ^ and a method of ^, we show that for any initial condition 
go in a measure class the (unique) solution q{t) of (II. ip is in a analytical functions space 
for all times t > 0. One consequence of the regularizing properties of equation (|l.ip is 
that many convergence phenomena of solutions towards equilibria or invariant manifolds 
happen in an analytic-functions space and not just in the classical L^ space. 

In section [3l we prove the existence of shadow inertial manifolds for the Kuramoto model 
(|l.ip . As inertial manifolds which have been introduced to overcome some fiaws of global 
attractors [12] [19], shadow inertial manifolds attract all solutions exponentially fast, the 
dynamics on a shadow inertial manifold is given by an ODE system, and each solution 
of the system has a phase on an asymptotically complete shadow inertial manifold. In 
that respect inertial manifolds and shadow inertial manifold give finite dimensional and 
accurate reduction of the long term behavior of a dynamical system. Our proof uses of 
original ideas of |27| . while avoiding some of its technical difficulties. 

These previous results establish essentially that the long term dynamics of equation 
(|l.ip can be captured by finite dimensional ODEs. From this point on, we focus on more 
specific properties of its asymptotic dynamics. In section HI we analyze the stability (or 
lack of stability) of equilibria. For readers' convenience, we start by briefiy recalling ex- 
isting results for equilibria of equation (II. ip of [8] S] mentioned above. Following this we 
perform local stability analyzes of equilibria. We write linearization theorems at the inco- 
herent equilibrium q{9) = 2^ in a rigorous mathematical setting, and we show that they 
confirm what has been largely expected in the literature. This thorough analysis also gives 
global stability information and estimates of escape times when the incoherent equilibrium 
is unstable, and it will be a basic ingredient in the study on the global attractor in section 
[H At the synchronized equilibria, an ad hoc Hilbert structure and associated spectral gap 
estimates have recently been found [^. Thanks to these we show that, when the coupling 
strength is large enough, K > Kc, the family of synchronized equilibria is asymptotically 
stable. Due to rotational invariance of (|l.ip . this family forms a circle like closed invariant 
curve. We give an estimate of the phase-shifting effect of a small perturbation of a synchro- 
nized equilibrium on this invariant curve. Finally, another important consequence of these 
local stability results is the uniform normal hyperbolicity of the family of synchronized 
equilibria that is of interest on its own. More on this can be found in |14) . 

In section[5l we establish several properties of the global attractor of (jl.ip . The existence 
of a global attractor in a classical consequence of regularizing properties of equation (jl.ip , 
such as proved in |31 [3 [5] . This is briefly recalled in section 15.11 In section 15. 2| relying on 
(un) stability analysis of equilibria in sectional we show that, just as without the diffusion 
term, equation (jl.ip possesses a two-dimensional invariant manifold. Interestingly, this 
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invariant manifold coincides with the global attractor of the model, whereas in the case 
of equation (ll.ip without diffusion, the global attractor, if it exists, cannot be reduced 
to the two dimensional invariant manifold of |21| : there are many unstable equilibria in 
the no-diffusion case that the global attractor would have to contain. Global attractors 
are classical tools to obtain a finite dimensional reduction of the asymptotic behavior of a 
large system [15j . this means that Kuramoto model (II. ip can asymptotically be reduced to 
a simple two-dimensional dynamical system on a disk. We characterize the trajectories 
on this manifold as orbits connecting the incoherent state to a synchronized one and 
further provide an algorithm for their numerical computation. Several properties of the 
global attractor of (II. ip . the dynamics on it and in its neighborhoods are also discussed 
numerically in section 15.31 In the strong-coupling limit (or equivalently in the small- 
diffusion limit ) K — )• -|-oo, we show that the global attractor of (jl.ip converges formally 
to the two-dimensional invariant manifold found by Ott and Antonsen J2T] in the no- 
diffusion case of equation (II. ip and supported by numerical evidences, we conjecture that 
this convergence holds in analytical functions spaces. 

In this way, we derive a full description at different levels of the long term dynamics 
of the mean field coupled Kuramoto model composed of infinitely many identical phase 
oscillators in the presence of noise. 

2. Existence, uniqueness and regularity results 

In the following we consider Sobolev spaces i/'^(S) = {g G ^^(S) : T.t'^ {^ + k^Y i^l + 1^1) < oo} 
for s > 0, and the Gevrey spaces g^ = [q£ L'^{S) : \\q\\l^ = Y.k=i "''("I + ^k) < +^} ■ 
When a > 1, the Gevrey space Qa is a subspace on the space of real analytic functions, 
for any k the natural injections Qa -^ H are continuous and compact, and if 1 < oi < 02 
then the injection Qa2 -^ Qai is continuous and compact. Notice that in a different context 
Gevrey functions are usually between C°° regularity and analytic regularity, but for partial 
differential equations the Gevrey spaces Qa above are classical too |11| . 

liq{t, e) = ^ + ^ Tl=i Xn{t) cos{ne)+yn{t) sin(n(9) is asolution of dLH) in C^{[0, -Foo[, L'^{S)), 
taking 2„ = x„ -|- iy„ one finds that for all n > 1 and t > 0, 

z'n = IZ-Zn + —[ziZn-l - 'zlZn+l] , (2.1) 

where zo{t) = 1 for all times. This ODE system is the key point to show that solutions 
are in a Gevrey space for any positive time, details of computations are similar to those of 
theorem 3.1 in [7] for the symmetric case y„ = (Vn > and t > 0). In [4J, it was shown 
that for Qo in a measure class, the solution q{t) of (II. ip is C°°(]0, -|-oo[xS"'^). Combining 
these, we can show that Gevrey regularization occurs in the general case (y„ 7^ 0) too, and 
regardless of the coupling strength K the solution of Kuramoto equation is in Qa for any 
a > 1 (for all t > 0). This proof also inspires the following lemma [Z!2| establishing that 
equation (|l.ip defines a continuous semiflow L'^ ^ Qa, which will be used several times in 
the next sections. 

Theorem 2.1. Let qo be a probability measure on S. There is a unique weak solution of 
the Kuramoto equation qt G C''([0, -|-oo[, A^(S)), absolutely continuous with respect to the 
Lebesgue measure for any t > 0, and denoting q{t, •) its density we have q{t, 9) > for all 
t > and 9 G [0, 27r]; we have q G C°°(]0, -|-oo[x [0, 27r]) and q{t,-) G Qa for any constant 
a > 1 and t > 0. Furthermore, for any a > 1 there is a bounded absorbing set Ba in Gevrey 
space Qa for equation M.l\) . 
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From now on we will denote St the semiflow associated to equation (|1.1|) . A bounded 
subset Ba is an absorbing set for a semiflow St, defined in a metric space X, if for any 
bounded set B G X the trajectories initiated in B enter Ba in finite time and remain in 
that set thereafter, ie 3to Vxq ^ B \/t > t^, x{t) G Ba- In theorem 12.11 above, Ba is a 
compact absorbing set for S't in X = H'^ for any s S N and in X = Qa' for any 1 < a' < a. 

The following lemma will be used several times in the following sections to show con- 
vergence results in Gevrey spaces. Its proof can be found in the Appendix. 

Lemma 2.2. (Gevrey- convergence Lemma) For any positive time e > and any constant 
a > 1, the semiflow S^ '■ L? ^ Ga associated to equation hl.l\) is Lipschitz-continuous 
on bounded sets. In particular, if q and q are two solutions of equation il-l\) bounded in 
L"^ , then for any e > 0, there is a Lipschitz constant L (depending only on e, Hgollia and 
\\qo\\L^) such that for any time t >0, 

Mt + e) - qit + e)|b, < L||(?(t) - mh^- (2-2) 

3. Asymptotically complete shadow inertial manifolds 

Inertial manifolds are finite dimensional, Lipschitz and positively invariant manifolds 
that exponentially attract solutions of evolution equations. There is an extensive litera- 
ture on the existence of inertial manifolds and related invariant sets such as inertial sets 
for partial differential equations ( |12| [23] [22] [T0| [9j). However, these general results do 
not apply directly to nonlinear Fokker-Planck equations similar to the Kuramoto model. 
Indeed, existence theorems for inertial manifolds require a spectral gap condition on the 
spectrum of the linearized equation that FPEs such as the Kuramoto model do not satisfy. 
This hurdle notwithstanding, Vukadinovic made a significant progress by establishing the 
existence of inertial manifolds for a number of non linear FPEs including the Kuramoto 
model ( |27| [28l 129)). For the latter, his result is stated in the case of even solutions. Our 
aim in this section is to show that a structure reminiscent of inertial manifolds exists for 
the Kuramoto model whether solutions are even or not. As some of the key steps of our 
proof are similar to Vukadinovic's work, we briefiy review his strategy. This short review 
also clarifies why we modify his approach in order to analyze the situation where initial 
data are not necessarily even functions. 

For a linear Fokker-Planck equation of the form 

dtq = dlq + de{qde{V)), (3.1) 

a standard method to obtain a regular reaction term is to set u = e^q, which gives 

dtu = dju + F[V]u. (3.2) 

Kuramoto equation (II. ip can be written in the form of (|3.ip . but is not linear because 
the "potential" V{9) = —Kj^ cos{9 — ip)q(ip)d(p depends on q. Nevertheless, in |27| 
Vukadinovic has shown that in the invariant subspace of even solutions of (II. ip . the map 
q I—)- e^q is a bi-Lipschitz bijection onto its image, and so that equation (j3.2p can be written 
in closed form in u. The reaction term u i— )• F[V{u)]u is then well defined H ^ H for the 
appropriate Hilbert space H, and standard theorems for the existence of inertial manifolds 
hold for (13. 2p . Using the inverse transform u = e^q i— )• q, Vukadinovic was able to establish 
the first inertial manifold existence results for equations of the form (|3.ip . Proving that 
the transform q \-^ u = e q is one to one, establishing its regularity and the regularity of 
its inverse are the most difficult and specific parts of |27) . 
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This method cannot be readily extended to the non-symmetric case, because it is not 
clear whether the map q i— )• V{q) is one-to-one for non-even q G L^. This motivates us 
to modify the approach to circumvent this question by embedding the system in a larger 
one where it is possible to prove the existence of inertial manifolds. This is sketched below 
with proofs left to the appendices. Prior to that we introduce the following definition of 
asymptotically complete shadow inertial manifold (AcSIM) and comment on it. 

Definition. Suppose that there exists a compact absorbing set Ba in X for the seniiflow 
St- Let Ai be a (finite Hausdorff dimensional) subset of X. The set Ai an asymptotically 
complete shadow inertial manifold (AcSIM) for St on X if 



• 



• 



there exists a smooth function <I> : M" — t- X such that Ai = <I>(1R"') is the graph of 

A^ attracts exponentially all trajectories of St-' there is a 5 > 

d{StXo,M) = o(e~^*) VxoGX, (3.3) 



• there is a flow St on M" (associated to an ODE system) such that and for any 
xq G X, there is a unique phase (pQ G W^ such that d{StXo, ^{Stcpo)) = O (e~ ) . 

An AcSIM gives a finite dimensional reduction of the dynamics of St^ to which the 
trajectories converge exponentially fast, and from a practical point of view this reduction 
to M gives as much information as an inertial manifold would do. The flow St is a shadow 
inertial form. An AcSIM is an (asymptotically complete) inertial set as soon as it is 
positively invariant by St- An AcSIM is a C^ manifold as soon as the function $ is C^ and 
injective, and if it is both a manifold and an inertial set, an AcSIM is a classical AcIM. 
The uniqueness of the phase (J)q implies that the trajectories PSt on M capture the slow 
dynamics of St on X (dynamics at speed less than O (e~ )). When they exist, inertial 
manifolds, inertial sets or approximate inertial manifolds contain the global attractor of 
a dynamical system. The same type of result holds for asymptotically complete shadow 
inertial manifolds. 

We introduce now several notations and we sketch the main steps of proof of the existence 
of AcSIMs for the Kuramoto model. Naturally, several of these steps are similar to [27] . 
notably when it comes to checking that the transformed system admits an inertial manifold. 
We consider the injection / : g i— )• (g, xi,yi) where xi = J^ q{9)cos{9)d9 and yi = 
/_ Qi(^) sm{6)d6, and we denote the projection P : (q, a, b) i— )• q. We define the transform 

J^ : a I— )■ W = a where V = —^ {acos{6) + bsm{6)). We will also denote 

\bj \ b ) 

u = e q- We will show that if g is a solution of equation (ll.ip . then U = J-oI{^q) is solution 
of a reaction-diffusion problem 

'^ + AU = N{U), (3.4) 

at 

where N is well defined H^ x M^ — t- H^ x R-^ for any s > 0. We will then check that equation 

(j3.4p is well posed, that it has a bounded absorbing set, and that the eigenvalues of A satisfy 

the spectral gap condition ensuring the existence of inertial manifolds M* for equation 

(j3.4p . We will denote by St the semifiow generated by equation (j3.4p in the following. 

Then, we will show that the transform T is regular and bijective, and we will deduce the 

existence of an AcSIM for equation (II. ip from the existence of M* . More precisely, we will 

show the existence of a finite dimensional C^ graph M. which attracts exponentially any 
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solution of equation (|1.1|) . For any initial condition go £ H^, the corresponding solution 
StQo of (II. ip has (for any large enough to) a unique phase ^o £ A^*, characterized by 

WSt+toQo - PJ'-'StVoWH^ = 0(e-''*), (3.5) 

where St is the restriction of S^ to (the finite dimensional manifold) M* . This graph Ai 
is an asymptotically complete shadow inertial manifold. 

This strategy makes use of a new kind of attractors, i.e. AcSIM. It does not rely on 
the injectivity of the potential q h^ V, which makes it simpler. It improves the functions 
spaces in which exponential convergence is proved from L^ to iJ^ for any s > 0, and it 
allows the definition of inertial form St = St\M* (the ODE system ruling the dynamics on 
Ai*), which allow numerical finite-dimensional approximation of the dynamics of (II. ip in 
principle. We will also see that the transform q i— )■ (e g, xi,yi) is C : H^ — t- H^ x M^ so 
that results of Rosa and Temam [22] may be applied to equation (j3.4p , and show that the 
inertial manifolds A4* (and the dynamics on it) are C'^ and uniformly normally hyperbolic. 
The proof of this theorem is detailed in the appendix. 

Theorem 3.1. For any s G N, the set Ai = PJ-^^ (Ai*) is a C^ finite dimensional graph 
in H" , there is a r] > and there is a flow St = 5'*|x* such that for any go £ H'^ and to 
large enough, there is a unique vq G A4* such that 

\\St+toqo - PJ'-^StVoWHs = O (e-^*) . (3.6) 

We have in particular distH'^{St+toQOi-M) ^ WSt+toQo — PJ^^^StVoWn'^ = O (e~^*) . The 
graph Ai is exponentially attractive for the flow St generated by equation U.l\) . StVo is the 
phase associated to the solution Stqo of IJA^, and PF~^StVQ is its shadow on A^. 

Furthermore, the set of equilibria of equation ( (g.^l ) is the image of the set of equilibria 
of (EIP by the transform T o I , and the global attractor of equation ^3.4^ is the image of 
the global attractor of (CUP by the transform T o F 

4. Equilibria of the system : Stability analysis 

4.1. On the stationary solutions to (II. ip . Up to a rotation all stationary solutions of 
(|l.ip which are probability densities - we call them equilibria - can be written as 

q{0) = ce^Kr cos(e) ^ (4 -^) 

with r a solution of the equation 

r = i^y^Ji , with L(s) = — [ (cos(e)y eMscos(9))d6 for s G M , (4.2) 

h{2Kr) 27r Jq 

and c such that j^^ q{6)d0 = 1. We refer to [8l|27] for this, and to [4j for this and a review 
of the mathematical physics literature on this issue, and a number of side facts. See |18| 
for similar results in the case of equation (11. ip with J * q{t, 9) = f^ sin {9 — ip)q{ip)d(p, 
and see [6j for asymptotic estimates when K — )• +oo in a three dimensional case similar to 

dnj. 

Equations (|4.ip and (j4.2p directly imply that the fixed point equation (|4.2p can have at 
most three solutions: r = 0, which is always present, and gives g(-) = 2^: this corresponds 
to the incoherent state of the system; and two non trivial solutions ±r (we denote by r the 
positive solution) that exist if and only if when K > 1. If we set qip{9) := q{9 + <p) , where 
g is given in (14. ip with r the nontrivial positive solution of (14. 2p . then g^ is an equilibrium 
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for every choice of ^j, so that we have a circle C = {q^ : (p £ S} of equihbria where each q,^ 
describes the (partially) synchronized state around the phase cp. 

Equation (ll.ip can be written in a gradient-flow form dtq = dg I qdg c fm I , where the 



>5q{e)J 

functional J^ is not increasing along solutions of equation (jl.ip . The existence of a Lya- 
punov functional J-" gives information about the asymptotic behavior of the system, and 
the stability of equilibria, see for example [2l [13] . 

In the following sections we will go beyond these results, by exploiting properties of the 
system linearized at the equilibria and by using invariant manifold theorems. Additionally, 
the following stability and unstability results will also be useful for the study of the global 
attractor (section [S]) . We do this separately for the incoherent stationary solution and for 
the synchronized solutions. 

4.2. The incoherent equilibrium ^. When K < 1, the constant ^ is the only sta- 
tionary solution, and the existence of a Lyapunov functional implies that it is globally 
attractive. 

Proposition 4.1. For any value of K and any initial condition go £ A^(S)) the co -limit set 
^(Qo) = 1^ {Stqojt ^ T} is included in the set of equilibria of equation U.l]) . In particular 

if K <1, for any qo G A^(S) we have \\q{t) — 2^|L — > 0. 

Proof. The convergence to g^p is a consequence of both the existence of a Lyapunov func- 
tional, and the fact that ^ is the unique equilibrium of (II. ip . Lemma [2.21 implies that this 
convergence happens in Gevrey spaces Qa- D 

We now use invariant manifolds to describe precisely the dynamics around 2:^ and to 
have estimates of the speed of convergence to 2:^ . 

Proposition 4.2. When K < 1, there are constants e > and M < 00 such that if 
1 1 go — 27rllw^ — ^ then the solution of equation ([HiP satisfies 

q{t) < M qo e"^* for every t > 0. (4.3) 

Moreover we have 

q{t) = ^ + P{qo)e-'^' + e(t), (4.4) 

with \\e{t)\\Hi < C\\qo - ^Ibie"'^* for any ^ < 7 < 1 - K and C = C{-f) > 0, 
and and P is continuous from a neighborhood of 2:^ in H^ into E i-k with the property 

P{q) = n(g— 2^) + \\\q — 2^\\fji) 1 where 11 is the orthogonal projection onto E i-k with 
kernel Ej^. 

Proof. The linearized operator Lj_v = i^d'gV — 2:^ J' * f is diagonal in the classical Fourier 

2.77 

basis of L^(§), and self-adjoint with compact resolvent. Classical proof of local stable 
manifold existence (see [17J chapter 5 section 1 for example) can be adapted here, where 
there are two identical largest eigenvalues instead of one only. D 

When K = 1 a bifurcation occurs at 2;;:. While 2:^ still attracts every initial condition 
(cf Prop. 14. ip . solutions do not approach y' exponentially fast. Nevertheless there is an 
exponentially attractive (in Gevrey norm) two-dimensional central manifold, along which 
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a phase is defined: for any solution q{t) of equation (|l.ip . there is a solution on the central 
manifold q{t) - the phase - such that q{t) converges to q{t) exponentially fast. 

Proposition 4.3. Suppose K = I. Let k > 2, k G'N. There is a map ^ G C'^{0,E+) - 
-E+ equipped with the H^ norm - with ^(0) =0 and D^(0) = 0, and a neighborhood O of 
the origin in E k-i such that the manifold 

2 

Mq = {uq + ^{uo) : no GO} (4.5) 

has the following properties: 

(i) The central manifold Aio is locally invariant: if q{t) is a solution of equation il.l\) 

on [0,r] with q{0) e O n Mq and q{t) E O for all t G [0,r], then q{t) € Mq for 

allt £ [0,T]. 
(ii) A4o contains all solutions which stay close enough to i^ for all time t G M; if q{t) 

is a solution of equation il-l\) for t G M such that q{t) G O for all t G M, then 

g(0) G A^o and q{t) G A^o for all t £R. 
(iii) A^o is invariant under reflections and rotations: if q & Mq then 9 i— t- q{—9) and 

6 I—)- q{0 + if) (for any (p G [0, 2-k]) are in M.q. 
(iv) A4.Q is locally attractive: for any qo G 7W(S), there is gg G Mq and constants C , 

5 > such that 

\\Stqo-St%\\g^<Ce-''. (4.6) 



Proof. The results (i) to (iv) are consequences of Theorems 2.9, 3.13 and 3.22 in |16j 
and the Gevrey-convergence Lemma (Lemma 12. 2p . These theorems apply in the setting 
Lj_ : H^ — )• L^ is linear continuous and the non linear part R G C^{H^,H^) where H'^ 
are the classical Sobolev spaces. D 

When X > 1, the incoherent equilibrium at g^ is no longer stable. Classical local 
stable and unstable invariant manifold theorems apply here, and due to the specific form 
of equation (II. ip . we also have more precise results about invariant manifolds, particularly 
about the way solutions leave ^. 

Proposition 4.4. Suppose K > 1. The stable manifold at i^ is 

W = W (^ =^ + E+ = {qeL^ : q>0, f q{e)de = 1, / q{e)e'^d9 = 0} . (4.7) 

(i) For any go G W^ , we have \\q{t) — ^\\g^ < C\\qQ — ^Wl'^g 2~* for all t > 0. 
(ii) Set O = Bl2{^,6) n {J^q{9)de =''l}. If 5 < 1 -\ and if a solution of (EJP 

satisfies q{t) G O for all times t > 0, then g(0) G W^ . 
(iii) More precisely, if for a 5 g]0, 1 — l/-ftr[ we have q{t) G O for every t G [0,T] then, 

with qo = q{0), we have 

I qo{0)e''^de = dist{qo, W) < Se-^^i^^^-^^-^^^ , (4.8) 

or, conversely, if \ f qo{9)e d9\ > e for a value of e ^ (0,6), then there is a time 
t < Ts^e = - K{l-5)-l ' ^^^^ ^^"^ 9(^) ^ ^L2(^,'^)- 
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Moreover for any a g]0, min(2, (X — l)/2[, there are constants C and a neighborhood 
O of w- in H^ such that there is a two-dimensional local unstable manifold Wf^^ Lipschitz 
continuous in H^ with the properties 

(iv) 2^ S W^ and W^ has a tangent space at g;^ which is E i-k . 

(v) For any uq G W^, equation l\l-l]) has a solution for t G] — oo, 0] with q{0) = 2^ +""0 
such that 



'<*' - i 



<C 



'(»' - h 



e"* , (4.9) 



for all t < 0. 
(vi) And for any qq ^ O with J qQ(9)d9 = 1, if there is a solution q{t) of equation ( fj.j)) 

which is defined on ] — cxd,0] and satisfies q{t) — > i^^, then qo G W^. 
The fact that ^^^"(2;^) is two dimensional will in particular be used in section [5^ 

Proof. The elementary, but key point is to realize that the dynamics on W^{J-) is linear: 
it is just dtu = ^f^e?- So (i) is a direct consequence of standard results on the linear heat 
equation or of Lemma 12.21 

To prove that any solution staying in a neighborhood of ^ is in the set W^^, that is 
(ii), we consider the Fourier ODE system (12. ID equivalent to equation (ll.ip . If we set 
Zn = Pn^^^" G C, we have 

P'n = -%-Pn + ^Pl [Pn-1 COs(6'i + 6n-l - On) - Pn+1 COs(-6'i + 6n+l - On)] , 
PnO'n = ^Pl [Pn~l Sm.{Oi + On-l - On) - Pn+1 sin(-6'i + 9n+l - On)] W . 

In particular, p'^ = — 2P1 + -fPi [1 ~ P2cos{02 — 20i)] . Now, observe that (7 G O implies 
Pn ^ 1 for all n > 1 and p2{q) < <5. So, if q{t) G O for all t > Tq, we have p'^ > 
^{K{l — S) — l)pi on [To, +00 [. Since 6 < 1 — ^ and i^ > 1 we have pi{t) > pi{To) exp(ct), 
with c > 0, so if Pi{Tq) > we have that pi{t) grows arbitrarily large, which is absurd, so 
we must have Pi{Tq) = 0, that is (/(Jq) G W^. 

For what concerns (iii) we use once again that q £ O implies p2 < S. Thanks to the 
Hilbert structure, we also see that dist{q,W^) = pi = [xf + yf)^'"^ ■ Now, suppose that 
q{t) G B{l/2-K,5) for ah t G [0,r]: like before we obtain p\ > pi\ {K{1 -6)- 1), and 

-K:(1-c5)-1 , K{\-S)-l , 

pi{t) > pi{0)e 2 . Therefore pi{0) < pi(0)e 2 . For the second statement in 

(iii) it suffices to choose pi(0) > e and to suppose that q{t) G B{l/27r,6) for all t G [0,T]: 
this leads to a contradiction for T sufficiently large. 

Statements (iv) to (vi) are classical results on the existence of an unstable manifold in 
a neighborhood of 2:^ (see for example |23|). See section 5 for further global results on the 
unstable manifold W^i-^). D 

4.3. The non-trivial equilibria. We assume K > 1. Several results of [3] will be 
used in this section, which we briefly sum up here for readers' convenience. Recall (cf 
§ 14. ip the notation(s) q{0) = qoiO) = ^^j }2Kr) '^'^^^^"^^^^ ^°^ °'^^ °^ ^^^ non-trivial equi- 
libria of (II. ip when K > 1. We consider the linearized operator Lq at q LgV = 2^!^ ~ 



Kdg [(J * q).v -\- {J * v)q] , and the scalar product <^ u,v ^= Jg ^^dO defined for u and v 
such that there are a Z^, V G L'^{§>, \dO) with W = u and V = v, and /g ^dO = /g ^dO = 0. 
The corresponding Hilbert space is denoted H^,.. It is not difficult to see that H^,~ is 
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equivalent to H ^, that is the space with weight one (see details on this issue in J14l 
Section 2]). 

Theorem 4.5. (cf\^) The operator L^ is essentially self-adjoint in H7',.. Its spectrum is 
pure point and lies in ] — oo, Ai] U {0}, where Ai < and is a simple eigenvalue of L~ 
with eigenvector dgq. Furthermore, for u,v £ D{Lq), we have 

1 I' uv f 

< LfjU, V >=< u, L^v >= -- / ^-de + v.{J* u)d0 , (4.10) 

and for K fixed we have D{LJ ) ~ L^ ,~(S) ~ -^^(S), meaning that the scalar products 
<C L^u,v ^, < u,v >= Jg^dO and {u,v) = J^uvdO are equivalent on R{L^) = {v G 
L2(S), f^vde = < v,deq >= 0}. 

Of course Ai = XiiK) depends on K. In the following we will content ourselves with 
the fact that Xi{K) < for every K > 1, see [4J for a sharper explicit bound on Xi{K). 

The next theorem shows that the only long term effect of a small L'^-perturbation of 
any synchronized equilibria ^^ G C is a small rotational shift. 

Theorem 4.6. There is 6 > such that if qo G L^ with qo > and f^qQ{9)d9 = 1 and if 
there exists ^ £ § such that Wqo — q^pWiz < 5 then then, for q{t) the solution of equation 
(QTTP with q{0) = qQ, we have \\q{t) — qipW^^ < 6 for all t > 0. Moreover there exists a 
(foo G S such that for any < /3 < | Ai | 

hit) - q^+^J\g. = 0{e-^') , (4.11) 

and we have 

\^oo- V>o\= oiWqo-qxl,^^) as \\qo-q^\\L2 ^0. (4.12) 

From the quoted theorem of [4j of from 14. 6| we can deduce that C is an invariant and 
stable normally hyperbolic curve for (II. ip . (see definition in |3| or [23] for instance). 

Corollary 4.7. The smooth and invariant curve C is stable normally hyperbolic. More 
precisely, at each point q^ G C, we consider the splitting 

L' = T^^C®N;, (4.13) 

with T^ C = M.dgcjip is the kernel of the linearized operator Lq at q^, and N^ is its 
orthogonal with respect to the H^,^ scalar product (TV* is also the range of the operator 
L^ ). This splitting depends continuously on y? and there is a T > such that 

\\DS^ {q^).v\\ < \\v\\ for any v G Tg^C and n G Z, and (4-14) 

2 

ll^'S'"^(^v')-^ll ^ ^ll^ll /o^ «^y ^ ^ %' «"^ n>l. (4.15) 

The proof of the theorem is outlined here, details and the proof of the corollary are 
postponed to the appendix. We mention that using that the invariant curve C is in fact 
not only invariant but a composed of equilibria only, its normal hyperbolicity can also be 
derived from the spectral gap result of [2] directly, see |14| . 

Proof. Thanks to rotation-symmetry of equation (II. Ih . without loss of generality we can 
assume that -0 = 0, ie q^ = g and \\qQ — q\\i2 < 6. 
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Lemma 4.8. Consider u{t, 9) = q(t, 6) - q{e) and x^{9) = q^{9) - q{9) = q{9 + ^) - q{9). 
There is a neighborhood V of in Lp' and to smooth projections 



2 



V ^ ]-6,e[cM V ^ Rm)cL 

u H> ip{u) u H> yu, 

such that u = x^{u) + y„ for all u G V, and u + q is a solution of equation \l-l\j if and 
only if ip = <p{u) and y = yu cli"g solutions of 

^ = $((^,y) and ^ + A2y = gi^,y) , (4.16) 

where 

H^^ y) = -z — J. — < V. f{x^ + y)- f{x^) > , (4.17) 

aiv, y) = f{Xf, + y) - f{x^) - d^x^<^{ip, y), (4.18) 

are smooth, with D^{0, 0) = and Dg(0, 0) = 0, and A2 is the restriction of L^ on R{L^), 
ie A2 = LglRi^L.y 

Proof. Details of proof of this lemma can be found in the appendix. D 

From this we can deduce a local stability result for the y variable. Let us introduce the 
natural Hilbert basis of eigenvectors e^ G D{Lq) so that L^ = 5]]fc>o ^k *^ "i e^ ^ e^ , with 
Ao = and > Ai > A„ > A„+i — ;■ —00 for n > 1. We have A2 = X]fc>i ^k ^ •, e^ 3> e^ 
and e*^2 = ^^>-,^ e*-^*-- <C •, e^ > e^, so that 

\\A^e'^H\\l., < e(i-^)^i* max (| A^^^e-^^l^'^-l*) \\v\\j,-^ , (4.19) 

for any e G]0, 1[ and from this we directly infer 

l|e*^^HlD(^|) < %^e-(i-)l^^l*||^|U-. (4.20) 

for any f > and a > (we have used ||'y|||)(-^a'i for ||^2'^||^-i). 

Lemma 4.9. For every /? g]0, |Ai|[ there exist po > and M > such that, if \ipo\ + 
II?/o||l2 = gP ^ gPO; then for all t > 

Mt)\ + \\y{t)\\L2<p. (4.21) 

Moreover for all t >0 we have 

\\y{t)\y<2\\yo\\L2e-'^K (4.22) 

Proof. Details of the proof of this lemma can be found in the appendix. D 

Since $ is smooth, we have -^(t) < C||yo||e^^*, and ip{t) converges: \(p{t) — (poo\ = 

0(e-^*). So we have 

\\yit)\\L2 + Ht) -Voo\ = (e-'^*) , (4.23) 

which, combined with the Gevrey-convergence lemma, gives (|4.1ip . 

For (|4.12p we argue instead that \-^\ < 2-j{p)\\yo\\i2e~^^ and that \(poo\ is bounded 
by 27(p)||yo||L2i = odV'ol + llyoll), for some -/{p) satisfying 7(p) — > (see the proof of 

lemma [^^ in the appendix), and o(|c^o| + II2/0II) = o (||go ~ ^||) since u 1— )• {<f,y) is Lipschitz 
and its inverse is Lipschitz too. This completes the proof of Theorem 14.61 D 
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5. The global attractor 

This section is entirely devoted to the study of the global attractor of equation (ll.lh . It 
starts with general results and progressively moves to more refined descriptions of this set. 
The end of the section proposes conjectures supported by numerical investigations. 

5.1. Existence. The existence of a global attractor is a classical consequence of the exis- 
tence of a compact absorbing set, and regularizing properties of equation (II. ip imply its 
existence (see [H 13 |5] for example). We state it here : 

Theorem 5.1. There is a set A bounded in Gevrey space Qa, such that for any space E 
such that the injection Qa ^ E is compact, for any hounded B set in E, we have 

dist{StB,A) — > (5.1) 

where dist is measured in the natural E norm. This includes the cases E = L?' , E = H^ 
for any s > 0, E = C for any k > I, and even E = Ga' if a' < a. 

When K < 1, the global attractor is {2;;^}. In the remainder of section [U we discuss the 
global attractor in the case K > 1. 

5.2. Description of the global attractor. Equation (jl.ip is equivariant under C(2), 
and the subspace of even functions is invariant by the dynamics. However for q S i^(S) 
there is in general no a S [0, 27r[ such that q{- + a) is even. More precisely, for q = 
27r ~'~ 27r Y2^kCOs{k9) + yf^sm{k9), such an a exists if and only if the quantities . ^'^ ^ 

and ^^ ^ are independent of A;. So the whole dynamics of (II. ip in L^ are not captured 

by its restriction on the subspace of even solutions. Nevertheless, we show below that the 
dynamics on the global attractor A are characterized by its restriction on the subspace of 
even functions. In fact the global attractor is radial: it is composed of one even heteroclinic 
solution qh, with qh{— 00) = ^ and qh{+oo) = q, and its rotations qh{t, • + 9?). 



Proposition 5.2. The unstable manifold ^^"(2:^) is the global attractor of the system. 

The 2 dimensional unstable manifold of the constant equilibrium ^^ consists in a family 
of heteroclinic orbits, each one connecting 1^ to a non trivial equilibrium q{- + ^p) (^ & 
[0,27r[j. All these orbits are obtained by rotation from one heteroclinic connection in the 
even functions subspace. In particular, for any solution q in W'^{2:^) of equation hl.l]) . there 
is an angle ip independent of the time t, such that q{t, ■ + ip) is even for any t G] — 00, +00 [. 

Proof. In dissipative systems unstable manifolds of equilibria are included in the global 
attractor [T5]. The converse inclusion is a rather general property, it relies essentially on 
the existence of a Lyapunov functional and on the fact that ^ is the only equilibrium with 
unstable directions. It implies that the global attractor A is exactly the unstable manifold 
W{^) here. 

We show now that in the even subspace, the 1-dimensional unstable manifold of the 
constant equilibrium ^ consists of two heteroclinic orbits, connecting 2:^ to q and ^(• + vr). 
Up to the symmetry 6 >-^ —6 and time translations, with no loss of generality we can 
consider only one solution q{t) escaping from 27 with xi{t) > when t — )• —00. The 
dissipativity and regularity properties imply that the trajectory {q{t), t G M} is compact 
and that q{t) it converges to one of the 3 equilibria of the system q, -^ or g(- + vr). One 
has the estimate |xfc(f)| < 1 for all Fourier coefhcients along solutions of our equation, and 
we can easily see that x\{t) > for all time f G R. From the Fourier ODEs system we see 
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that lim xi{t) = is absurd if xi is not for all times. Hence the solution q{t) is not 

t— >+oo 

an homoclinic connection of ^, and it is a heteroclinic solution. The 0{2) invariance of 
equation (jl.ip eventually implies the second part of proposition 15.21 D 

5.3. Global attractor and numerical explorations. As in the previous section, we 
denote q^ the heteroclinic solution of (jl.ip such that qh{t) — > ,- and qh{t) — > q, 

i— s>— oo t— !>+oo 

where q{9) is the unique equilibrium of (jl.ip with xi > in the even functions subspace. 
Since the structure of the global attractor is radial, without loss of generality, we can 
restrict ourselves to working in the subspace of even functions. 

This section is organized as a sketch of a proof of two conjectures, which are partially 
proved and partially supported by numerical evidences. The first one is that the global 
attractor is not only a two dimensional manifold homeomorphic to a disk, but it is an 
analytical graph over the disk {xf + yf < 1} in L^, and an recursive procedure for the 
Taylor series coefficient is given. The second conjecture is that this graph, which depends 
on the coupling constant K, converges in analytic functions space when K — )• +oo to 
Ott-Antonsen two-dimensional invariant manifold. 

The global attractor is an analytic graph. Consider the global attractor as a geometrical 
curve T-L = {(?h(t), t €] — oo, +oo[}. 

Proposition 5.3. There is a neighborhood V of -^ in L^ such that TiCiV is a graph on 
Spanjcos}. That is to say there is e > such that 

nnv ^ ]o,e[ 

q ^ ^^ = J^^q(^g^cosie)de ^^■^> 

is a smooth bijection. 

Moreover, if t i— )• xi{t) = J^ qfi{t,6)cos{9)d9 is increasing on ] — cxd,+oo[ (qf^ denotes 
the heteroclinic solution of U.l\) ). then Ti is globally a C°° graph on Spanjcos}. 

Proof. By unstable manifold theorem at f^, we know that 7^ is a C°°-graph over M. cos 
in a neighborhood of ^, and there is a Tq > — oo such that t i— )• xi{t) is increasing on 
] — oo,To[. The first part of the proposition holds in particular for V = {q £ L'^, < 
J- q{^) cos{6)d6 < xi(Tq)}. If we have Tq = +c«, then 1 1— )■ xi{t) is a bijection from M to 
its image, and Ti is globally a graph on Spanjcos}. D 

In the following we denote the graph function by <I>-^ • "i inverse of the 



bijection (|5.2p . 

Numerical result. For all t ^M., x'i{t) > 0. (see figure [K 



Xl I— )• 



The property x'„(t) > is numerically true for any n > 1 and any time t € M along 
the solution qh. All these curves have a unique inflexion point, they are convex for large 
negative times t and concave when t — )• +oo. Their limit, the Fourier coefficients of q, are 
given by the Bessel modified functions : lim x„(t) = ^(2x1) ' ^° ^^^^^ ^^^ °f them tend to 

lim -f4~Y = 1 when K — )■ +oo, but for any finite K we have lim |xri(i)| < C4r for some 
constants C > 0, a > 1 and all n large enough. 
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Figure 1. The heteroclinic solutions q^ is a graph on xi. Plots of xi{t) 
along the heteroclinic solution q^ for K — 2.0 (bottom) and K — 4.0 (top). These 
curves were obtained using the Fourier ODEs system equivalent to (|l.ip , truncated 
with N — 200 equations, Euler scheme and dt = 0.00001. Numerics were checked 
using N — 100 and dt = 0.0001 and classical Runge Kutta forth order scheme, the 
maximum difference between the two is of order lO^'"^. The first Fourier coefficient 
xi is increasing along the heteroclinic solution. The same holds for all the Xn{t), 
n > 1 (Figures not shown). 



Recall that J_ q{9) cos{9)d9 = r (see section (|4.ip ). so that lim J_ qh{0) cos{t, 
liin xi{t) = r. If proved, the monotonicity of xi would imply the following. 

t— >+oo 

Conjecture. The graph function <l>-^ is C°°{[0,r]), and the global attractor W^ is a C°° 
graph on the disk {(xi,2/i), x\ + yl < r}. 

With the notations of proposition 15.31 above, there are functions (/?„ £ C°° such that 
for all xi G]0,e[, we have Xn{t) = J^^ qh(t, 9) cos{n9)d9 = ipn{xi{t)). So for t £] — oo,ro[ 
equation (|1.1|) is equivalent to 



x[ 



1 K 

--xi + —xi{l-ip2{xi)) and x^ = c^„(xi), Vn > 2. 



(5.3) 



Proposition 5.4. The functions ipn are C°° in a neighborhood of 0, their Taylor series 
are of the form 

ipn{x) = x" ^ an,2px'^^, 
p>0 
where the coefficients a„^2p can be computed recursively. 

Proof. We use the notation (pi{x) = x, so that ai^ = 1 and ai^p = for all p > 1. We 
use the Fourier ODE system x^ = — ^x„ + ■^xi(x„_i — Xn+i) (Vn > 1), equivalent to 
equation (II. Ih . and the Taylor expansion of the (/?„ to obtain relations between the a„^ 



p- 
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First one finds that p < n implies ifn = 0, so that (/?„,(x) = x" ^ an^pX^^, and then all odd 
coefficients an,2p+i are zero, so that the series has the form ipn{x) = x"' '^„>Q0Cn,2px'^^ ■ 

Then one finds a2,o = 'kTi a-^id 0^,0 = j^ 1 fa„i) Q^n-i,o = n"'~2 ttT? ^^^ ^^^ recurrence 
relations 

((1 + p)K + (1 - p))a2,2p = y Yl (2 + 2j>2,2ia2,2i - i^a3,2(p-i) (5.4) 

yi+j=p-i,i>o,j>o y 

for n = 2 and all p > 1 and 

((n + 2p)i^ + (n^ - n - 2p)) a„,2p = 

(5.5) 




n 



2i)an,2ia2,2i + Kn{an-i,2p - an+l,2(p-l)) 



for all n > 3 and p>l. The no-resonance conditions (n + 2p) + " ~^~ ^ 7^ are satisfied 
as soon as i^T > 1. 

We introduce the order ^ on N^ : we have {mi,qi) <C {rn2,q2) iff {qi < Q2 or (gi = 52 
and mi < m2))- This is the lexicographical order for the inverted couples {qi,mi). One 
can check that for that order, the recursion relation defining a„.2p uses only am,2q with 
{m,2q) <^ {n,2p). There are no infinite strictly decreasing sequences in N^ for the order 
<^, and so one can solve recursively the relations defining the an,2p. D 

Consider q^ an heteroclinic solution of (jl.ip and Xn{t) its Fourier coefficients. Since 
t I—)- xi(t) is increasing and bijective R — 7']0,r[, we can define Xn{x) for all x G]0,r[ by 
Xn{xi{t)) = x„ (x5"^(xi(t))) for ah t £ R. 

Numerical results. For any K and any n, the Taylor series defining ipn has a positive 
radius, which is r (see figure \^. For any K and any n, the sum of the Taylor series 
defining ipn equals Xn on ]0, r[ (see figure\3(j. 

For a series ^^ OpX^, the convergence radius R is given by R^^ = limsup|ap|p. In 

p— >'+00 

1 
figure [21 we see that r„„ = log((a„.2p)p+i) — > I > 0, where / > is a real number 

p— >+oo 

depending on K but not on n. That is to say all series ^ an,2px'^'^ have a positive radius 

R = e~''^, which depends on K but does not depend on n. We also observe that for 
all K, the limit is / w — 21og(r), suggesting that the radius R = r. Since the Taylor 
series of (/?„ are x" ^ an,2px'^'^ the same results hold for them. We have also checked that 
d'Alembert's ratio test is numerically consistent with this: for all n > 2 and iT > 1 we 
have ""'^^^^ — y I. and / ~ -^ suggesting that R = r too (figures not shown). Figure [3] 

shows that the troncated sum of the Taylor series (^„ ~ ^q an,2p2;"'~''^^ equals X„, on ]0, r[. 

Notice that the verification Ylp'^n,2pX^~^'^^ = Xn{x) is necessary since some functions, like 

_ 1 
X I— )• e ^ , have converging Taylor series without being analytical. 

The previous numerical observations suggest the following. The graph function satisfies 
'^nix) = — + -"Y ^Pnix) cos(n-) = ^ + - X] p" X] "ri,2px'^^ cos(n-), 

n>l n>l y p=l J 

and it is analytical on [0, r[ in the sense that all ipn are. 
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Figure 2. Positive convergence radius of the series X]p>o '^n,2pX^^~^"^ ■ Plots 



of: 



n.p 



— ]-r log(a„,2p) for different n and different values of K 



(a) : for n = 2, 



and from top to bottom K = f .5, K ^ 2.0 K ^ 3.0, (b) : for i^ = 2.0, and from 
top to bottom a = 5, a = 2a — 4. Notice that (a) and (b) show Vn^p in different 
windows. The radii of the series are positive since the sequences r„.p converge to 
positive values, with limp_j.+oo ?'n,p — ~21og(r). For any value of K, the radius 
R{ipn) does not seem to depend on n. The radius increases when K increases, but 
the convergence Vn^p — > — 21og(r) is slower for large K. 



Figure 3. 




sum 



E 



the function 

-1, 



^p •-^ii,,2pX 

Xn = Xn{xi). Plots of functions y — (p2{x) and y{x) ~ X2{xi^{x)). The 
sum of the series was computed with N — 200 terms, and the drawn curve 
{(a;,2/(x)), X € [0, r]} — {(xi(t),X2(t)), t G M} was computed with the Fourier 
ODEs system with TV = 200 equations, Euler scheme and dt — 0.00001. The nu- 
merical difference between the two is less than 10~^ for xi G [0,r]. The equilibria 
q IS, dX X ~ r^ characterized hy f2{f) = ^ — -h = \ here. The sum of the Taylor 
series coincides with the function Xn — Xn{xi) along the heteroclinic solution qh 
for all n > 2 . In particular, this heteroclinic solution seen has a graph on xi is 
an analytic function on [0, r]. 
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Conjecture. The global attractorli is an analytical graph on the disk {{xi,yi), xf + yf < 
r}. 

These Taylor series can also be used to estimate the eigenvalues at the incoherent equi- 
librium 2^^: and at the synchronized equilibrium q. At 2:^ one can analytically compute the 
eigenvalues, and we have checked that the numerical estimates given by the Taylor series 
are accurate. At q, the eigenvalues can also be computed by linearizing the truncated 
Fourier ODE system (12. ip . and we have checked that the two methods agree there too. 

The recurrences equations (15. 4p and (j5.5p can actually be used to compute recursively 
the Taylor coefficients a„^2p, and so one can compute directly the heteroclinic solution 

Qh{xi) = 7— + -y^^(pn{xi)cos{ne) 

and the whole global attractor with arbitrary precision without discretizing the partial 
differential equation (jl.ip or the equivalent infinite Fourier system (|2.ip . 

Convergence to Ott-Antonsen Ansatz for large K. 

Proposition 5.5. Considering the Taylor series (pn{x) = x^ "^ ^q an,2px'^^ , we have 
Q^n,0 — ^ 1 C'lT'd Oin,2p — > for all Ti > 2 and all p > 1. 

A'— >-+oo ' A'— 5>+oo 

Proof. The formulas a„,o = n"~2 'k+j found above show that a^.o — ^ 1 when K — )• +00. 
When p = 1 we have 

// 2\ n^ — n — 2\ 

( ( 1 + - j H — j an,2 = an,oa2fl + {an-1,2 - an+1,0) (5.6) 

and knowing an,o -^ 1, the recursion hypothesis an,2 — ^ for all n < A^ implies aAr+1^2 — ^ 0. 
One can directly check that 02 2 — ^ 0, and so deduce an 2 — > for all n > 2. 

When p > 2 (and n > 2), using what we know about the case p = and p = 1, and the 
relation 

^^^ 1 + 2-1 an,2ja2,2i 




Z-^ \ ^ 



•a.- 1 V -/ / (5.7) 

'^n+l,2(p-l)) ) 

we see that the recursion hypothesis 

y{m,2q) <C {n,2p) with q> I, am,2q — > 

implies a„ 2v — ^ and so this is true for all p > 2 and n > 2 by induction principle. D 

In the no-diffusion case of equation (jl.ip . the graph graph TioA = {QOAix), x € [0, ![} C 
L^, where qoA{x) = 27 ^ tF Si °° x" cos(n0) (x G [0, 1[), is an invariant manifold for the 
dynamics (see plj). In that case, the two-dimensional invariant manifold of Ott and 
Antonsen is obtained from qoA by rotations 

1 1^°° 
A^OA = {7^ + -yx"cos(n(6' + a)), xe [0,1[, aG [0, 27r]} C L^ 
zvr vr ^-^ 

1 

Proposition 15.51 implies that (fnix) — > x" and that <I>-^ — t- qoA formally. 
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Proposition 5.6. Assume that for all K large enough, we have x'^{t) > for all t G 
] — oo,+oo[. Consider P/v : L^ — ^ Span{cos{k9) , k < N} the natural L'^ -orthogonal 
projection. For all z g]0, 1[ and N < +oo, we have 

Pn ^niz) — > Pn QOAiz) with convergence in Qa for any a > 0. 

Proof. Recall that (see section [4.ip r = [^ q{9)cos{9)d9 = lim xi{t) — > 1. For 

any K large enough, we have xi : M. — )-]0,r[ is a bijection. For any a G]0, 1[ and K 
large enough, we have < a < r, and there is a unique ta^K ^] ~ oo,+oo[ such that 
JZ-^Qh{ia,K-,&)cos{9)d9 = xiita^K) = CL. The previous proposition implies that 

P^$^(z) = - + -j; Yl «n,2p^"+'M cos(n-) ^^ - + -^."cos(n.), 

n=l yo<n+2p<Af / n=l 

uniformly on [0, 27r]. The degree oi Pn ^h{z) (trigonometric polynomial) is independent 
of K, and so convergence occurs in any Qa, a> Q. D 

We recall now a classical property of analytic functions. 

Proposition 5.7. Let fx (K G Nj 6e a sequence of analytic functions with radius Rk > a 
for some a > independent of K . Suppose that g is analytic with radius Rg > a and 
fx — > g uniformly on [0,a]. Then we have fx^gin analytic function space, and 

in particular for any I > 0, -ih:fK — ^ -ih:9 uniformly on [0,a]. 

As explained above, figures [2] and [3] suggest that c^„ are analytical for any n, and their 
Taylor series X^j, an,2pa;^^ have the same radius R = r — > 1. 

Numerical result. For any n we have ipn{x) — > a;" uniformly on any compact subset 
o/[0, 1[ (see figure]^. 

Figure U shows convergence (pn — t- x" uniformly on any compact of [0, 1[ when K — )■ +oo. 

The difference (/?„ — Xn is increasing on [0,r] (for any K and n), and at x = r we have 

lim (pnif) — r" = 0, which shows that the convergence is uniform. Furthermore we 

expect theoretically \ipn{r) — r" = 0| = -^ + o(-j^) when X — )■ 0, and this is confirmed 
numerically on figure H] too, suggesting that the convergence happens with speed i. 
Assuming that these numerical observations are true, we have 

1 1 ^ ( 
^w(^) = ^ + - IZ p'' Yl «n,2p2^^ I cos(n- 

n=l y 0<n+2p<+oo 

with X^o<n+2p<+oo ^n,2pz'^^ — > 1 Uniformly on any compact subset of [0, 1[ (for any n), 
and I X^o<n+2p<+oo '^n,2pz'^^\ < 1 for all n, K and z G [0, r]. This would imply |99„(z)| < z" 
(for all n > 1 and z G [0,r]), and ^niz) = -^ + - X]n=i fn{z)cos{n-) is in Gevrey space 
Qa uniformly on K, for all z < a and a < —. 

Conjecture. For any a G [0, 1[ and any z G [0, a], we have 

^h{z) — y qoA{z) m Qa Vl<a<-. 
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Figure 4. Uniform convergence c/3„(x) — > x". Plots of functions i/)^ (a;) 
for X e [0, r], for different values of i^ : K = 1.5 (a), K = 2.0 (b), K = 4.0 (c), and 
K = 16.0 (d). All values were computed with N = 100 terms in the entire series 
defining (p2 ■ The top most curve in all plots is y = x^ . The functions (p2 converges 
to the limit y = x'^ when K — > +cx). uniformly on any closed set of [0, 1[. The 
difference x" — '^n{x) is increasing for x € [0, r[, so that supp^^.^^ jx" — ipn{x)\ = 
r" - fn{r) = (/ q{e) cos{e)de)" - J q{0) cos{ne)de. One has r^^l ^ ii? + O {^) 
and (finif) = 2 I (2K ) ~ ^ ~ Tk ^ ^ (k^)" Thus one expects a convergence 
towards the Ott and Antonsen two dimensional invariant manifold at speed -^r : 



^UPo<2;<r F 



fn{x)\ < ^ 



The two dimensional global attractor ^^"(2;^) of il-1]) converges to the two dimensional 
invariant manifold AioA of Kuramoto equation with no diffusion, in analytic functions 
space, with speed -^ when K — )• +00. 

The dynamics on T-L are given by x'^ = —\xi + Kxi{l — tp2{x\). Considering xi{t) = 
Xi{j^), in the limit K — )• +00, we find x'l = xi(l — xf), which is the equation ruling the 
dynamics on Mqa is the no diffusion case. 

Conjecture. Appropriately time-rescaled by a factor -^, the solutions of equation il-1]) 
on W^{2^) converges to solutions of the no-diffusion case on A4oa, in Gevrey space on all 
intervals ] — 00, T], with speed -^. 
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Appendix 

Proof of lemma 12.21 We will use the classical fact that for any positive time r/ > 0, the 
semiflow Sn of equation (jl.ip is well defined L-^ : — )■ L^ and Lipschitz continuous on bounded 
sets, and the Lipschitz constant can be chosen uniformly for r] in bounded intervals [0, e] 
(see |17) for example). 

Consider q^ and (fo in L^, and q{t) and q{t) the associated solutions of equation (jl.ip . 
The equation for the difference w = q — q is linear : 

dtw = ^d^w - Kde{Fw) - Kde{{J * w)q) (5.8) 

where F{t) = J* q{t) = xi{t) cos(6') - yi{t) sin(6'). 

Noting w = ^Yl ianCos{ne) + /3„ sin(n6')), 7n = "n + if^n, zi = xi + iyi and Zn = 
Xn + Wm the previous equation implies the Fourier system : 

In = -y^n + —{Zi-fn-l - zTjn+l) + — (7l^n-l -Tl^n+l) (5-9) 

with 70 = for aU t > 0. And taking An = zi^'^n-i + ^Jnj^ and B„ = -^ij;;zn-i + 
yijnZn-i we have (similarly to the previous proof) 

^|7nP = -n^\ln\ + — (^n " ^n+l) + — (^n " ^n+l)- (5.10) 

Notice that for any time r/ > and for any n > 1, 

\ln{ri)\ < 1^(77)11^2 = \\q{ri) - q{r])\\L2 < Lip{S^)\\qo - qoh^, (5.11) 

so for /(n) = a"™™(*.i)^ using |zi| < 2, we have 

df^fin), ,o\ ( 21og(a) C a^-A/J^ ^^^, .A , , 

M 

<Cf{N)\^N\\lN-i\+C\-ii\ Y. /(n)|7„||z„_i| (5.13) 

n=Af+l 



<C(t)||c 



9o\\h + N^N + l) (£ ^f(^)\^nA+C\j^\' f ^ n/(n)|5„|2j (5.14) 



C ^ '' 



< C{t)\\qo-qo\\l2 + j^jY—r I Yl ^/(^)l7n|^ ]+C{t)\\qo-qo\\l2\\q\\l^Mi,t)^N,M (^-l^) 



\n=N 
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From theorem 12. H we have that |Q'(t)|G' ^i^/n) is bounded on [0,e], for any a > 1, and 
so ||Q^(t)||G minfi t) ^^ bounded on [0,€] too (for any a > 1). The constants C above are all 
independent of M, so there is a iV large enough such that for all t S [0, e], 

1/2 / \ 1/2 




fi'^h.^ |2 1 , ^ J y^ ^^^^M.. |2 



7nr I +C\}_^ n/(n)|7„|^ I < C{t)\\q,{0) - 92(0)^2. (5.16) 

n>N 



This proves that for any qq in L^, for any e > and 5 > 0, for any a > 1, there is a 
M < +00 such that for any qq G 5^2(^0,6) we have \\q{5)-q{5)\\G^^^„^^s) ^ ^lko-^o||L2- 

Proof of theorem 13.11 Definition of the transform T and the new equation (P^. In the 



following we slightly abuse notations by using the same letter V to denote both 

V =V{q){9) = -Kixi{q)cos{e)+yiiq)sm{9)), 
or V =V{a,b){d) = -K{acos{e) + bsm{0)) ^ ' 

depending on context. For q{t) a solution of (II. ip . the ordinary differential equations 
satisfied by xi = xi{q{t)) and yi = yi{q{t)) are x'^ = —■^xi + -^ [xi{l — X2) — ^1^2] and 
y[ = "2^1 + Y ["-^1^2 + 2/1(1 + X2)] and this gives in particular 

dtV = --^ - f [Fi(n,xi,yi)cos(^) + Gi{u,xuyi)sm{e)] . (5.18) 

Along a solution of the Kuramoto equation (|l.ip . where u = e q and 

Fi{u, a, b) = [a(l — X2{u, a, b)) — bY2{u, a, b)] 
Gi{u,a,b) = [-aY2{u,a,b) + 6(1 + X2{u,a,b))] , 

with 

X2{u,a,b) = j-^^^,(6/)ef («™«W+''«i>^W) cos(20)d^, 
Y2{u, a, b) = Z^ u{^oy^{acos(e)+bs\n{e)) s[n{29)d9. 
We have that for any solution q of equation (II. ID , U satisfies (13. 4p where 



(5.19) 
(5.20) 




AU= \ 1 XI = |xi , (5.21 



K ri 
' 2 




[IV + 1(99^)2 + f Fi(^/) cos(0) + f Gi(Z^) sin 



iV(Z^) = ( hKF^iU) I (5.22) 

iKGiiU) 



We will show now that inertial manifolds exist for equation (13. 4p . There are several 
inertial manifolds existence theorem (see |26 1 \29 \ [22] for instance), such as the following. 

Theorem 5.8. Consider an abstract evolution equation on a Hilbert space H 

dii 

— + Au = f{u) (5.23) 

where A is a positive self-adjoint operator with compact resolvent and f : H ^ H is Lips- 
chitz continuous on bounded sets. The spectrum of A consists in a sequence of eigenvalues 
Xn such that A„ < A„-|_i — > +00. Suppose that equation Ii5.23\) has a bounded absorbing 

set in H . Then there exists a constant C , such that if 

A„+i -\n>C (5.24) 
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for some n > 0, then there exists an asymptotically complete inertial manifold Ai* for 
equation Ii5.23\) . which is the graph of a Lipschitz function ^ on En = (Bi<n ker(j4 — Xild). 
Furthermore, if f is C^{H,H), then the graph function $ and the inertial manifold Ai* 
are C^ too. 

During the next two steps of this proof, we show that theorem 15.81 apphes to equation 



The Cauchy problem for equation l^cl.4^ is well-posed. The operator A is diagonal in the 
natural Hilbert basis of E := H^ x M x M, and the spectral gap condition of theorem 15.81 
holds for the operator A. The maps X2,Y2,Fi, and Gi : L^ x M x M — > M are C°°. Since 
the first line of the non linear term N{U) is the product of a Fourier polynomial (whose 
coefficients are C°^{H' x M x M)) and u, we have that TV is C~ : H' xR"^ ^ H" xR"^. In 
particular, the non linear term A^ is Lipschitz on bounded subsets of any H'^ x M x M, and 
classical semi-linear parabolic equations results in Hilbert spaces apply (see for example 
0). 

Proposition 5.9. Let s G N. For any Uq G H^'^'^ x M x M, there is a unique solution U{t) 
of ll^P with U £ C°{[0,TmaxbH'+^ X M X M) n C^{[0,Tmax[; H' X M X M), and we have 
either Tmax = +oo or lim ||W(t)||H=xRxlR = +oo. 

(- ^ ^ "max 

Furthermore, ifUo > 0, then we have U{t) > for all t G ^,T„iax[- 

In the following we denote by S^ the semiflow generated by equation (j3.4p . Thanks to 
the regularizing properties of equation (|l.ip . without loss of generality we may assume that 
qo G H^ and Uq(^ H" xRxR (for any fixed s > 0). 



Existence of a hounded absorbing set for equation ^3.4\j . We have already established 
that all hypotheses of theorem 15.81 hold for equation (|3.4p . but the existence of a bounded 
absorbing set. 

For any lA solution of (j3.4p . {q, a, h) = T^^ (U) satisfies 

dtq = \dlq + Kde [(-asin(-) +hcos{-))q] , (5.25) 

a' = -ia + ii^Fi(e-^('*=°'^(-)+''^iO)g,a,6), , . 

b' = -i6 + iKGi(e-^(«-«(-)+^«''^(-))g,a,6). ^^^^> 

As for Kuramoto equation (jl.ip , we prove that equations (|5.25p and (j5.26p have a bounded 
absorbing set in a Gevrey space, by mean of an equivalent Fourier ODE system. 

Denoting wi{t) = a{t) + ib{t) G C, and z„(t) = x„(t) + iy„(t) = f^^e^"-%{t,0)d0 we 
have 

n Kn 

z'n = --;rZn + —[wiZn-i -wTzn+i]- (5.27) 

With the method of theorem 12.11 and lemma 12. 2| we obtain 

(5.28) 
If wi is bounded on some time interval [0, T], then for A^ large enough we have ^1^1 f at + 
Ikll/AT < C, and for any solution U G C^{[0,T],L'^ x M x R), we have u{t) G Hf 
for all t g]0;T]. We also see from the previous argument that if we have a uniform 
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bound |u;i| < Ci for all t G [0, +oo[, then N can be chosen large enough such that 
^klf AT + Ikllf Af — C* and the bounded ball i?// (0, 2C) in Gevrey space is an absorbing 
set for our system. 

We now prove a uniform bound for wi on [0, +oo[. Equation (|5.25p preserves the integral 
of q, J^ q{t)dO = J^ qodO along solutions. Since go ^ implies q{t) > for solutions of 
equations (j5.25p and (j5.26p . we have |xi| < J_ qo{0)dO = C and |yi| < C as long as the 
solution exists (where xi = J_ q{t,6) cos{6)d6 and yi = J_ q{t,6)sm{9)d9). Consider 
{qo,ao,bo) G L^ x M x M and the associated solution {q(t),a{t),b{t)) of equations (j5.25p 
and ([5T2B|) . We have 



x[ = -ixi + i/^Fi(e-^("=°^(-)+''^^'^(-))g,a,6), 



Vi 



' - iyi + iKGi(e-^("^°"(-)+*^^'^(-))g,a,6), 



(5.29) 



so that "i~(^i ~ '^) = —[xi — a) and — (yi — b) = —{yi — b). (5.30) 

(JjL (JjL 

The quantities \xi — a\ and |yi — 6| are exponentially decreasing. Since |xi(t)| < 1 and 
|2/i(OI — 1 fo'^ ^11 time i > 0, there is a time to such that \wi\ < \a\ + |6| < 4 for all t > to- 
As seen in the previous paragraph, this implies that 

^kl/,iv + MIn <Ci Vt G [0, +oo[, (5.31) 

and the bounded ball i?// (0,2Ci) in Gevrey space is an absorbing set for our equations 
(f5l25]) and (fOHjl and so for equation (|3^ . 

Existence of an AcIM for equation ^3.4^ and an AcSIM for equation lil.l\) . By theorem 
15.81 there exists an inertial manifold M* for equation p.4p . Furthermore, since the non- 
linear term J\f {U) is C'^ (A: > 1), the inertial manifold M* is C^, and for any Uq G H^ x M? 
and to large enough there is a unique phase Vo G A^* such that 

\\S:+,^^Uo-S:Vo\\=0{e~^') (5.32) 

and the leaves 

£vo = {^o, \\S:Uo - S^VoW = O {e-^^')} (5-33) 

form a C^ foliation of space H^ x M?. 

To deduce properties of the flow of equation (jl.ip from results for solutions of equation 
(|3.4p . it is essential that the transform F o I : q \-^ U and the inverse J-~^ are smooth. 
With the definition u{9) = g((9)e^^^(^i'2'i'^), one sees that 

II^IIl^ < Ci^lkllLH|e-^'^("^^°^(-)+^^^'"")||L- < CKC{K,xi,yi\\qh2) (5.34) 

and also for any s G N ||ii||_ffs < Cx,sC{K,s,xi,yi)\\q\\H'', so that g i— t- n is well defined 
H'^ — )• H^. Besides the mappings q i— )• (xi,yi) i— )• V{xi,yi,-) i— )• e~ are smooth and 
one finally finds that F : q i-^ U is C'^{H^;H^ x M x M) for any s G N. For the same 

reasons, the inverse mapping F^^ U = {u,a,b) i— )■ {q,a,b) = (uea^^"'^"'^^'-*^'"^™^'-*^) is also 
well defined and C~ smooth F'* x M x M ^ i7* x M x M for any s > 0. 

For any solution q of equation (jl.ip . U = F o I{q) is a solution of (|3.4p . and {(g, a, b) G 
H^ X R^, xi{q) = a, yi{q) = 6} is an invariant manifold for equation (13. 4p . To any go ^ H^ , 
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we associate Uq = J-iqo), which has a phase on the inertial manifold Ai* of equations (|5.25|) 
and (I5.26p . ie there is a ig ^ind a Vq S A^Jj such that 

\\S*{t + to)Uo - S*{t)Vo\\e = O (e-"*) . 

Consider the projection P : H^ x M'^ — >■ H'^ defined by P{q,a,b) = q. Taking v{t) = 
PJ-^^S^Vq, since the projection P and the inverse J^~^ are smooth, we have 

\\S{t + to)qo-v{t)\\Hs=0{e-'^'). 

The inertial manifold Ai* is finite dimensional, so that Sf = S*\_m* is a fiow on a finite 
dimensional space. Since Ad* is a smooth graph on E^ (see theorem IS.Sp . At = PJ-~^Ai* 
is a smooth graph too. We conclude that At = PF~^Ai* is an AcSIM for equation (jl.ip . 

Proof of lemma 14. 8L To simplify the notation we will write H~^ as as a short-cut for 
H^,. in this proof. Recall that Lq is normal, and we have the orthogonal sum L? = 

L>(L|)=ker(L^)ei?(L^). 

Let us denote by P the orthonormal projection onto ker(L|j). For u and x^p in a neigh- 
borhoods of the origins, the equality u = Xip + y with y G R{Lq) is equivalent to Pu = Px^ 
and y = Zu — z^. We check that for each u in a neighborhood of there is a unique x^,(„) 
such that Pu = Pxi^r^\, and variables change u i— )■ (x^(„), Zu — 2<^(u)) will be well defined. 

Since x^ = q{- + c/?) -q{-), and -^|^^q = dgq / 0, the set {x^ : ip g] -e,e[} is a graph 
above R • d^q = ker(L^) for small enough e. In particular that for u close enough to zero, 
there is a unique ^{u) such that Pu = Px^(„). Moreover, (/9 i— )• x^ is C (] — e,e[,H^,~) 
{k > 0) and the decomposition u i— )• (a;^(„),y = ^i^ — -z,^(n)) is C''^ in a neighborhood of 
u = in H^, . The mapping u i— t- (c/?, y) is well defined and C : H^, — t-] — e, e[x/i'~, on 
a neighborhood of the origin. 

We can now write the dynamics for the new variables cp and y. The x^ are equilibria, 
ie we have L^x^p = f(x^), for ip G [0, 27r]. If u is a solution of -^ + L^u = f{u), we have 

(d^Xp) -^ + -^ + ^qV = fi^f + y)- fi^f) ■ (5-35) 

Let us denote A2 the restriction of L^ on R{L^), ie A2 = L^jl^/i^^, and take v = doq/Wd^qW'^ 
(so that L^v = and <C v,dipXip ^ |(^=o =^ v,dgq ^= 1). Then we get the system 

^ = $(99,2/) and ^ + A2y = giip,y), (5.36) 

where 

^('^' 2/) = :p — P, — ^ ^ ^' -^(^^ + ^) " -^(^^^ ^ ' (^■^'^) 

and 

'/(V', y) = fixip + y)- f{x^) - d^x^^{ip, y) . (5.38) 

We notice that the non linearity / is well defined i^(S) — )• H^^{E>), and since / is 
bilinear, we have / G C^{L'^,H~^) (like before, k is arbitrary). Since / G C^{L'^,H^^) , 
Xip G C ([0, 27r]; L^) and 3<^a;y,|<^=o = deq 7^ 0, it is easy to check that there is a e > and 
an ?7 > such that $ :] - e, e[xB£,2(0, 77) — )■ M and g ■] - e, e[xi?j;^2(0, r/) — > if"-*^ are C'^ . 
The Frechet derivative of / at the origin is null, ie Df{0) = 0, and the expressions for $ 
and g give D^{0, 0) = and Dg{0, 0) = 0. 
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Proof of lemma [4.91 Inequality (I4.20p also holds if we put -D(^2 ") ^^^ -^(^2) instead 
of 15(^2) and H^, . For a = 0, one directly has ||e^^u||£)(^7) < e"''* 11^11^/^7-) for any 
t > 0. We now choose a = 1/2 and a value for e. 

Since the Frechet derivative of $ and g are smooth, for any 6 > 0, there is an 77 > 
such that \(p\ + ||y|| < r] implies 

\\D^^,y)\\<5, and \\Dg{ip,y)\\ < 6 , (5.39) 

where the norms are the appropriate operator norms. We also see that for any 99 s] — e, e[, 
we have '&((/7, 0) = and g{(p, 0) = 0. Thus for any /? > 0, there is a j{p) such that on the 
set {((/7,y); \^\ + WvWl^ < p}, we have 

\'^iv^,y)\ + Mv, y)\\H-^ <i{p)\\y\\L^, and 7(^)^0. (5.40) 

By the hypothesis on the initial condition and by the continuity of 1 1— )• I'/'tl + ||yt||_L2 we 
have that there exists T > such that \ipt\ + ||yt||L2 < p. We check first that the estimates 
of the lemma are true on [0,T]. 

For yo £ Ri^q), we have 

yit) = e-*^2yo + / e-^'-'^^^gi^is),y{s))ds , (5.41) 

Jo 

which implies on [0, T] 





t 



(5.42) 



< l|yo||L^e-l^^l* + Ci/2,e / ^=^e-(i-)l^^l(*-^)7(p)||y(5)||L^ds, 

Jo Vt- s 

where e is arbitrarily chosen in ]0, l/2[. 

Then we consider u{t) := supo<^<( ||y(s)||^2e'"'^^^'^'"^'*. We have for t G [0, T] 

||y(i)||i2e(i-2^)l^il* < llyoIlL^ +7(p)C7i/2,e r^=^e-^l^^l(*-^)n(s)d5 

Jo Vt-s 

<\\yo\y+lip)MMT), (5.43) 
where M^ := Ci/2,e /n°° ^e~'^'^'^*^ds and this means (choose p so that 7(/9) < l/(2Me)) 

n(r)<2||yo||L2, (5.44) 

which directly implies (^^?m with /3 := (1 - 2e)|Ai|, still for t G [0,T]. 

We now want to relax the condition t G [0,r]. For this we first observe that since 

1^1 = \'^{^,y)\ < l{p)\\y\\L2, we have 

dip 
'dt 

and therefore \(p{t)\ < Iv'ol + 7(p)^- If we now choose po such that we have jip) < /3, 
then on [0, T] we have 

\ipit)\<\ipo\ + \p<lp, (5.46) 

and 

\\y{t)\\L2 + \^{t)\ < 2\\yo\\L2 + Ip<Ip<P- (5.47) 



<2j{p)\\yo\\L2e-^', (5.45) 
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This implies that there is a Ti > T such that ||y(t)||L2 + \f{t)\ < p on [0,Ti] and therefore 
that T can be chosen arbitrarily large. In fact if we set T^, := sup{i : |</?t| + lll/tllL^ ^ p} E^nd 
if Tj. < oo we can take T = T^ in the first part of the proof and we arrive to a contradiction. 

Proof of corollary 14.71 Because of rotation symmetry, it is enough to study the system 
around q = g^=o- We consider the splitting E = Lp' = ker(L^)©i?(L^) = T^CQNL setting 
N? = R[Lq). This splitting extends continuously to all q[- + 93) G C by rotations. 

Any point of C is an equilibrium, so that for any time t and integer n, we have Snt\c = 
Id\c, and so DSnt{S[)-v = v for any v G T^C. This implies that ||L)S"-^(g).r;|| < \\v\\ for any 
T > and n G Z. 

We consider now (9701^0) ^■iid {(p{t),y(t)) the associated solution. The estimates from 
the lemma 149] above imply \ip{t) — ipo\ < 27(p)4||yo||i and ||y(i)|| < 2||yo||e~^*, where 

p = 8(|(^ol + l|yoll)- 

For V £ R{Lq) = N~, we have v = (0, y^,) = (0, v) in the new coordinates, and if we 
consider q{0) = (930, yo) = q + ev, we find 

'Mp) 



Ht)\ + \\ym< 



/3 



+ 2e 



-lit 



\yo\ 



(5.48) 



With llyoll — sP ~ ^ll^ll' '^^ have ^{p) — > 0, and so there is a ei such that for any 



e < ei we have ^)f' < -j^. And choosing T such that e " < Tjn gives 



S' {q + ev) - S' 



< 



1 




1 1 




+ 2 




iU 




2uJ 



20 



e\\v\\ < -\\v\ 
II II - 5II I 



(5.49) 



and \\DS^{q).v\\ < i||w||. 

In the same way, for any n > 1, there is a e^ > such that for all e < e„ < ei, we have 



g •^ < ™r, and then 



^nT 



TIT- 



S'"' {q + ev) - S"' q\\ < 



1 1 

+ 2 — 

10*^ 20" 

'nTf 



e K" < \\v\ 

II II - -^Q^ll I 



(5.50) 



for any v G R{Lq) = N^, and this implies HDS"" (g).!)]! < yqttII^II- This completes the 
proof of normal hyperbolicity of C. 
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